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THEORY  OF  THE  TRANSVERSE  GRADIENT  HIGGLER 


I.  Introduction 

J.  Madey  and  co-workers  have  proposed  the  use  of  an  FEL  with  a 
transverse  field  gradient  to  prevent  the  energy  spreading  and  detuning 
which  saturates  the  ordinary  FEL.  The  basic  idea  is  that  the  transverse 
field  gradient  acts  as  a  spectrometer  so  that  as  an  electron's  energy 
changes  due  to  Interaction  with  radiation  It  is  focused  to  a  different 
transverse  position.  One  then  arranges  the  transverse  field  gradient  so 
that  the  resonance  condition  k  ■  ^  k  remains  satisfied  at  all 

•  777^  » 

tranverse  positions,  x  .  Here  k8  is  the  radiation  wave  number,  kw  che 

eR. 

nolav  MOWA  numkAW  on/)  a  _  **  f  Ka  aw  nnf  4  a  1  a  f  a  **4  aaI  aw 


wiggler  wave  number,  and  a 


,  the  vector  potential  of  the  wiggler 


magnet  field,  i8  usually  close  to  unity  in  value.  We  will  concentrate 

our  attention  on  the  possibility  of  obtaining  high  extraction  from  such  a 
magnet. 


Madeyl  has  proposed  a  wiggler  structure  characterized  in  the  y  »  0 


plane  by: 


S(z,x)  -  ey  |V2B(x)  cos  kwz  +  Bc(x)j 


B(x)  -  B, 


[(1  +i)2 


2  +  x/t 
.  2 


and 


Bc(x) 


1  e  B(x)B'(x) 

^omcV  1+i 
w  n 


Thus  the  field  gradients  are  characterized  by  n~*  .  This  choice 

preserves  the  transverse  Independence  of  resonance.  The  equations  of 
motion  for  an  electron  orbiting  In  the  x  -  z  plane  [after  averaging  over 
the  fast  (kw“*)  motion]  have  been  given  by  Madey  as: 
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x 
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(3) 


Here  |E|  Is  the  electric  field  of  the  optical  wave,  6  -  X-  -  1  wittl 

Yo 

YQ  being  the  design  Initial  energy  for  a  particle  which  enters  at  x  ■  0, 
and  ^  is  the  optical  phase  ^  M  •flt  .  q  ,  the 

optical  phase  slip.  Is,  of  course,  a  function  of  frequency.  The  design  has 
been  chosen  such  that  v^  ,  and  hence  q  is  Independent  of  electron 
energy  so  long  as  the  betatron  focussing  keeps  the  electron  at  x  ■  * 

Further,  we  have 


2 


and  J  is  a  constant  close  to  unity 


The  underlined  terms  in  equations  (1)  -  (3)  are  small.  Moreover, 
there  is  some  disagreement  among  Madey's  papers,  and  further  with  recent 
work  of  Morton,  as  to  their  exact  coefficients  which  depend  on  details  of 
design.  However,  all  workers  are  agreed  as  to  their  approximate  magnitude 
and  form  and  we  shall  see  that  the  nature  of  the  solution  is  insensitive  to 
the  details  of  these  small  non-linear  effects.  We  will  first  consider  the 
linearized  problem,  neglecting  the  underlined  terms. 

First,  however,  it  is  convenient  to  introduce  dimensionless 
variables.  Let 
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x>2Tnx 

v 


5  +TXT* 


q 


(5) 


The  equations  transform  to  (keeping  terms  of  order  —  ) 
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X 
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-(x  -  6) (1  -  —  S) 
k  - 
w 
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e  sin  ip 
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(8) 


It  may  be  seen  that  the  non-linear  terns  are  all  formally  of  order 

k  e2|E|B  a  a  k2  2 

«  1  •  Further  c  =  -  -  -  J  ■  4  — — —r  — ^  J  =  — =■  where  £1  is  the 

k  2  4,  22  ,,2,2  ,2 

w  me  k  y  1  +  a  k.  k. 

B  o  w  B  B 

synchrotron  frequency  for  oscillation  in  the  ponderomotive  potential 

wells.  In  Sections  II  and  IV  we  discuss  regimes  in  which 
2 

e  *  (fl/kg )  «  1  .  The  case  e  »  1  is  discussed  in  Section  V. 
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II.  Linearized  Solution  for  e  <<  1  . 

The  linearized  forms  of  equations  (6),  (7),  and  (8)  are: 

x  -  -(x  -  6)  (9) 

S  *  e  sin  $  (10) 

$-q-(x-d)«q  +  x  (II) 

Equation  (11)  may  be  integrated  immediately  to  yield 

$  *  ♦q  +  qz  +  *  (12) 

where  we  have  absorbed  the  initial  value  of  x  into  the  initial  (random) 
optical  phase  <Jiq  .  Differentiating  equation  (9)  then  yields  the  desired 
equation  for  X  =  x  . 


X  -  -X  +  e  sin  (i|>  +  qz  +•  X)  (13) 

o 

Equation  (13)  is  evidently  of  the  form  of  a  driven  harmonic 
oscillator.  We  are  interested  in  achieving  a  high  degree  of  excitation. 
Thus  the  interesting  regime  is  q  ■  1  .a  nearly  resonant  oscillator.  We 
look  for  a  solution  of  the  form 

X  ■  f(z)  sin  [qz  +  *(z)]  (14) 
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r 


where  £  and  $  are  slowly  varying  in  z  . 


Thus 


••  •  •  r  ••  •  2  1 

X  •  2(q  +  $)  f  cos  (qz  +  $)  +  j^f  -  (q  +  $)  f  I  sin  (qz  + 


♦  > 


-  2f  q  cos  (qz  +  $)  -  f(q2  +  2q  |)  sin  (qz  +  $)  (15) 

In  deriving  equation  (15)  we  have  made  use  of  the  assumption  that 
f  and  $  are  slowly  varying. 

Further,  using  the  Bessel  function  expansion 


If  sin  X  In  X  T  ... 

!  "  2-* e  vf) 


we  have 


sln(t(»o  +  qz  +  X)  -  Jn-1^  sin  [n^qz  +  <fr)  ”  (<fr  ~  i|»q)J 


Here  we  need  only  the  secular  term  (n  *  0)  and  the  resonant  terms 


n  *  #*• 


Thus 


6 


sin  (<|i  +  qz  +  X)  «  J.(f)  sin  ($  -  *  ) 


+  2J|  (f)  sin  (qz  +  $)  cos  ($  -  tj^) 


J^f) 


cos  (qz  +  $)  sin  ($  -  i|>q) 


We  are  now  In  a  position  to  substitute  equations  (IS)  and  (16) 
into  (13)  and  equate  coefficients  of  sin  (qz  +  $)  and  cos  (qz  +  $) 
First,  however,  we  specialize  to  the  near  resonant  case  by  defining 


1  -  q2  «  qX 


and  setting  q  =  1  elsewhere  in  the  equations.  From  equation  (13),  we 
then  have 


,  Ji<f) 

f  “  -e  — _ —  sin  ($  -  \J>  ) 


f|  mY-  -  eJ'(f)  cos  (*  -  *q) 


Further,  we  may  determine  the  energy  change  from  equation  (10)  and 
the  secular  term  in  equation  (16). 


6  ■  eJ,(f)  sin  )  *  -ff 

1  o 


so  that 


f2.  f2 

«  -  . 

2 

with  f  «  1  the  initial  value  of  f  . 
o 


(20) 


Note  that  the  non-secular  terms  in  equation  (16)  would  lead  to 
oscillatory  terms  in  6  which  are  0( c)  «  1  ,  indicating  small  energy 

spread. 


Dividing  equation  (19)  by  equation  (18)  we  have 


d  cos  ($  -  ^»Q) 

- 3! - 


Xf  J1 


with  solution 


cos  (♦  -  V  -  TJ-  [c  +  h  fZ]  •  (22> 

We  will  be  concerned  with  low  initial  emittances  f0  ^  *  ln 
which  case  C  -  0  .In  this  case  we  see  that  near 
f  -  0  ,  cos  ($  -  $q)  -  0  ,  while  from  equation  (19)  it  is  clear  the  the 

stable  point  is  ♦  ”  “  y  •  Hence  initially  sin  ($  -  ^Q)  •  -1  and 

we  see  from  equation  (18)  that  f  initially  increases.  We  may  now  write 
down  the  solution  to  equation  (18) 


8 


(23) 


Equations  (20)  and  (23)  represent  tha  solution  to  tha  problem.  Saturation 
will  occur  for  f  given  by 

8ftX 


J1  <faa*> 


Since  e  la  by  assuaptlon  small  we  find  for  particles  off 

2*  2 

resonance  (X  ♦  0)  that  f  - and  -6  -  —Z—  .  For  resonant 

“X  *  “X  / 

particles  and  an  infinitely  long  wiggler,  growth  will  stop  at  Jj(f)  ■  0  , 

fi  o\2 

l.e.,  f  -  3.8  .  Here  8  ■  -  v  1  '  .  Thus  the  maximum  extraction 

max  2 

which  can  be  obtained  by  such  a  gain-expanded  wiggler  In  the  small  c 
regime  Is  (using  equation  5): 

-  e-l  -  (3.8)2  k8  (24) 

T  ^  k. 

By  comparison,  the  maximum  extraction  from  a  simple  wiggler  Is 
Q/k^  ,  which  Is  smaller.  (The  saturation  length  for  a  simple  wiggler  Is 

5  •> 

We  may  obtain  an  explicit  evaluation  of  equation  (23)  in  the  small 
signal  limit  since  Jj(f )  -  -jf  for  f  small.  Then  from  (23)  and  (20)  we 

have 
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(25) 


«Y  kB  2  kg  c2  2  X 

"Y--zirf  " IT "T iin  7*  * 

v  v  X 


A  reasonable  approximation  to  the  departure  from  the  small  signal 
limit  for  X2  >  3c2  la  obtained  by  approximating  J*(f)  -  f^(l  -  f^)  , 
which  yields: 


fiY  _  kB  e2 


,  2 

.in2  t2  *  f  i  . 


(25' ) 


The  shape  of  the  resonant  (X  ■  0)  extraction  curve ,  obtained  by 
numerical  evaluation  of  equation  (23)  is  shown  in  Figure  I.  An  approximate 
saturation  length  can  be  defined  by 


8 


r 


FIGURE  1  Oscillation  amplitude  f  at  resonance,  shown  as  a  function  of  interaction  length 
for  the  intermediate  Intensity  regime. 
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III.  Non-Ideal  Effects 


t 


> 


» 


i 


t 


t 


» 
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We  iiy  now  discuss  Che  sensitivity  and  modification  of  our 
solution  by  various  perturbing  effects. 

(1)  Horizontal  Emlttance 

We  have  assumed  that  at  a  ■  0  ,  x  -  n6  and  x  »  0  .If 
this  is  not  the  case  then  the  constant  C  of  equation  (22)  may  no  longer 
be  set  to  zero.  A  complete  analysis  in  the  small  signal  limit  at  resonance 
is  straightforward  and  shows  that,  after  phase  averaging,  <S>  is 
unaffected  by  the  value  of  fQ  even  when  f^  >  f^  -  f^  ,  so  long  as 
fQ  <  1  .  It  is  also  clear,  however,  that  as  the  initial  fQ  increases 

significantly  beyond  unity,  the  saturated  limit  of  6  decreases.  We  may 
estimate  that  the  horizontal  emlttance  Eh  is  related  to  fQ  by 


so  that  the  condition  on  maximum  allowed  horizontal  emlttance  is 

2  'S*  - 

f  -  2E  k  _  <  1  (27) 

o  h  s  k 

6 

(2)  Vertical  Emlttance 

For  a  simple  2-D  gain  expander  of  the  type  proposed  by  Madey 
the  vertical  and  horizontal  oscillations  are  decoupled.  We  have  not 
investigated  the  possibility  of  a  more  complex  3-D  structure.  For  the 
decoupled  magnet  it  may  be  seen  that  the  vertical  oscillation  frequency 
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The  variation  in 


kgy  is  of  the  same  order  of  magnitude  as  kg  . 
flutter  motion  with  y  makes  a  comparable  contribution  to  Av  •  The 
exact  relations  depend  upon  details  of  magnet  design*  As  a  simple 
parametrization  we  write 


phase  slip  be  less  than  *  and 
find 


2  2Av 


v  k'  k' 

0  0 


where  k.  •**  k'  .  Requiring  that  the 
0  0 

using  equation  (26)  to  estimate  L  we 


(28) 


(3)  Diffraction 

The  lengths  implied  by  equation  (26)  may  also  pose 
difficulties  in  passing  a  diffraction  limited  beam  through  the  magnet* 
Using  the  criterion  that  the  optical  beam  Intensity  for  a  Gaussian  beam 
should  be  down  by  e^  at  the  magnet  aperture  from  its  central  value*  and 
using  equation  (26)  to  estimate  beam  length  we  find  a  restriction  on  beam 
radius 


r 


2 

b 


(29) 


Since  geometrical  considerations  may  also  require  r,  ~  n 

0 

equation  (29)  may  provide  considerable  restraint  on  magnet  design* 
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'  rvr.— -r  -  ^ 


Equations  (26)  -  (29)  and  the  definitions: 


k  -  —  k 

8  TT7  w 

w 


.  i  ~+ 4 


q2  _  4asaw  u2 

u  - - 2  w 

1  +  aw 


eBW 

a..  -  — «=- 


me 


(30) 

(31) 

(32) 

(33) 


provide  the  basic  lmput  for  specifying  the  parameters  of  a  medium- 
extraction  single-pass  gain-expanded  magnet. 
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IV.  Effect  of  the  Omitted  Non-Linear  Terms 

We  turn  now  to  a  consideration  of  the  modifications  Introduced  by 
the  underlined  terras  in  equations  (6)  -  (8).  First,  however,  we  briefly 
consider  a  related  problem,  the  effect  of  using  wiggler  magnets  whose 
properties  ,  k^  ,  etc.,  vary  with  z  .  It  seems  clear  that  the 

principal  effect,  in  view  of  the  sharp  resonance  behavior  exhibited  by 
equation  (25),  will  be  to  make  A(z)  variable.  Thus  if  we  modify  our 
basic  ansatz  (14)  to  assuming  a  solution  of  the  form 
X  »  f(z)  sin  (z  -  y  J \dz  +  <fr)  we  may  reduce  the  problem  to  a  solution  of 
equations  (18)  and  (19)  with  A  replaced  by  A(z)  .  For  general  A(z) 
these  may  be  reduced  to  a  non-linear  2nd  order  differential  equation  for 
f  .  However,  it  is  clear  that  the  saturation  limit  of  equation  (24)  will 
not  be  exceeded. 


Let  us  now  return  to  the  non-linear  problem  and  consider  equations 

2  2 

(6),  (7),  and  (8).  Working  to  order  x  ,6  ,  x6  ,  etc.,  and  making 

use  of  the  fact  that  5  is  slowly  varying  to  0(e)  while  x  has  a  large 
oscillatory  piece,  we  may  rewrite  equation  (8)  as 


■  q  +  x 


1  +  (x  -  5  +  26) 
w 


k8  *2 

4k  X 
w 


q  +  X 


1  +lf5 

w 


B  r  *2  .  "2  \ 
-  4ir  (x  +  *  ) 


w 


so  that 
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♦  ■  q*  +  ♦, 


-  qz  +  ^ 


■^/Ci 

-  *ir  +  51 


+  X  )dz  +  X 


1  +w~ 6 

w 


xZ  )dz  +  X 


(34) 


and  from  equadon  (6) 


X  “  -(1  6)  X  +  6  (1  -JL  6) 


(35) 


Again  we  make  Che  ansaCz 


X  ■  f  sin  ij>  +  qz  - 


[*o  +  qZ  "  ^ Ji*2  +  x2)dz  +  ♦] 


and  keep  only  those  terms  which  affect  the  resonance.  Thus  it  is  adequate 

k„ 


to  use  S  -  e  sin  <|»  ,  neglecting  the  order  modifications  in  e 

w 

Equations  (18)  and  (20)  are  then  unaltered  while  equation  (19)  becomes 
[recalling  that  &  -  f  sin  (  );  X  ■  f  cos  (  ) ,  6  -  -f2/2  ] 


fl  2 

x+r 

w 


eJJ(f) 

— - - cos  (#  -  <|>  ) 

f  o 


(36) 


As  we  have  done  in  the  linear  theory  we  may  integrate  equations 
(18)  and  (36)  to  obtain 


cos  (*  -  <|io) 


1  TX  f2  +  kS  f4l 


(37) 
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and 


if 


fdf 


JL  *2  -  6 


IT 

v 


/Be)’ 


1/2 


Xz  v  w 

Recalling  that  the  resonance  width  is  determined  by  —  <  —  , 

4  2 

kB  -  4. 

we  see  that  the  non-linear  effect  will  be  unimportant  when  _  <  _  . 

k  2 

w  k  Lf 

2  0 
Taking  f  »  8  as  a  typical  saturation  value  and  using  (26),  this  criterion 


,  K0  -  4it 

may  also  be  written  —  <  __  £  .  When  this  criterion  is  well  satisfied 

k  64 


w 

the  saturation  level  of  equation  (24)  becomes  accessible  if  X  is  given  an 
appropriately  chosen  small  negative  value. 


We  see  that  non-linear  effects  become  important  when  — is 

w 

large  enough.  For  the  resonant  electron,  X  ■  0  ,  peak  extraction  will 

then  occur  for 


or 


at 


A  k 

f_ _ 0 

8e  k 

w 


Jj(f) 


f/2 


(38) 


z 


r2(7/6) 
r  (4/3) 
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By  detuning,  i.e.,  choosing  a  slightly  negative  X  ,  a  factor  of 
about  2.S  improvement  can  be  obtained  in  extraction,  while  of  course  a 
carefully  chosen  variable  pitch  X(z)  could  cancel  the  non-linear 
frequency  shift  and  recover  the  saturated  efficiency  of  equation  (24).  As 
discussed  earlier,  different  magnet  designs  could  also  result  in  slightly 
different  non-linear  terms  thereby  yielding  a  different  numerical 
coefficient  in  equation  (38).  However,  it  is  very  difficult  to  see  how  a 
design  of  this  type  could  improve  on  the  extraction  efficiency  given  by 
equation  (24).  We  have  considered  in  this  note  only  designs  for  q  -  1  as 
proposed  for  small  e  high  single  pass  extraction.  However,  the 
techniques  employed  herein  may  also  prove  of  use  in  analyzing  the  behavior 
of  functionally  canceling  designs  proposed  for  storage  ring  adaption. 
Finally,  we  note  that  operation  at  q  ■  0  could  also  give  comparable 
extraction  if  a  pre-bunched  beam  (Klystron)  were  used. 

We  have  Integrated  equations  (6)  -  (8)  numerically  for  the  case 

k 

e  *  .05  ,  i— ■  .2  for  which  equation  (38)  predicts  (-  — )  »  .050 

Kw  y  max 

Numerically  we  obtained  (-  Al)  -  .051  . 

y  max 


V.  The  High  Intensity  Regime 

We  recall  from  the  theory  of  the  uniform  wiggler  that  the 
parameter  ft  corresponds  to  the  small  amplitude  oscillation  frequency  in 
the  ponderamotive  potential  wells,  and  may  be  referred  to  as  the  optical 
synchrotron  frequency*  In  our  previous  discussion  we  have  assumed  this 
frequency  to  be  small  compared  to  ,  the  betatron  frequency,  here 

induced  by  the  transverse  gradient*  The  case  in  which  ft  ~  k  can  be 

M 

expected  to  be  complicated  as  it  is  likely  to  lead  to  coupling  between  the 
two  oscillation  modes.  However,  the  case  ft  »  k.  again  leads  to  behavior 

P 

which  can  be  discussed  in  a  simple  way. 

We  first  differentiate  equation  (11)  to  obtain  t|i  ■  -(x  -  6)  , 

and  using  equations  (10)  and  (12)  we  find 


■  e  sin  <J>  +  -  <|r  +  qz 


(39) 


To  relate  the  discussion  to  that  given  in  Ref.  (2)  we  shift  the  phase  in 
(39)  of  both  ip  and  i|»o  by  w  .  That  is,  we  write  ^  «  $  +  n  f  which 
yields 


$  ■  -e  sin  $  +  +  qz  (40) 


n2 

Since  e  "  — ^  ,  we  are  here  Interested  in  the  region  e  »  i  .  For 

qz  «  c  ,  §he  equation  is  approximated  by  ^  »  -e  sin  f  .  Since 

|  (0)  -  q  (for  x  ■  5)  ,  those  electrons  with  -*  +  ^  <  *  -  ■!=. 

y/c  /e 
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are  crapped  in  oscillatory  orbits*  That  is,  all  but  a  fraction  q/« 
are  trapped,  and  we  assume  this  quantity  to  be  very  small.  For  the  trapped 
orbits  1^1  ,  |$q|  <  v  and  we  can  therefore  consistently  treat  these  terms 
as  small  compared  to  e  •  The  term  qz  ,  of  course,  grows  with  z  and 
we  shall  show  below  that  it  leads  to  a  progressive  detrapping  of  particles, 
but  it  also  leads  to  a  reduction  of  particle  energy  for  those  which  remain 
trapped.  Appealing  to  the  theory  of  the  variable  parameter  wiggler 
discussed  in  Ref.  2,  we  shall  assume  that  the  main  contribution  to  the 
averaged  energy  exchange  comes  from  the  trapped  particles,  and  neglect  any 
contribution  from  particles  which  have  become  detrapped.  Thus  for  trapped 
particles  equation  (40)  may  be  written 

^  £  sin  ^  +  qz  *  ^  (41) 

where 

F($)  ■  -e  (cos  ♦  +  ♦  sin  ifr^)  (42) 

with 

sin  (43) 

Equation  (41)  describes  motion  in  the  potential  F($)  ,  which  has  the 

general  form  of  Figure  2.1  of  Ref.  2.  The  trapped  particles  execute 
oscillatory  motion  about  the  minimum,  which  occurs  at  ^  a  •  and  for 
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rrzzi 


the  trapped  particles  <sln  $>  ,  which  represents  the  average  over  this 

motion.  Is  equal  to  sin  .  This  statement  is  rigorously  true  only  for 

z  independent  ,  but  the  same  result  holds  approximately  when  we  may 

regard  as  varying  adlabatically .  Again  assuming  adiabatic  variation, 

the  phase  area  of  any  trapped  orbit  remains  constant,  but  the  maximum  phase 
area  available  for  trapped  orbits  shrinks  as  Increases  and  vanishes 

as  $r  reaches  j  .  There  is  thus  a  continual  detrapping  of  particles 
as  z  increases,  and  when  qz  ■  e  all  particles  are  detrapped.  Since 
there  is  then  no  further  average  energy  exchange,  the  wlggler  should  be 
terminated  when,  or  before,  this  point  is  reached. 

On  the  basis  of  the  above  physical  description  we  write 

<6>  »  -e  sin  4/  jCsin  i)>  ) 
r  r 

-  -qz  j(-|^)  (44) 

where  j(sln  $r)  is  the  area  In  phase  space  enclosed  by  the  last  trapped 
orbit  relative  to  that  at  $  ■  f)  .It  is  given  explicitly  by 

*  a(sin  *  ^-)  defined  by  equation  2.60  in  Ref.  (2).  It  follows 
that  j(o)  is  unity,  and  if  ^2.  <<  1  ,  equation  (44)  yields  the  simple 

result 


1  2 
<«>  "  -  2  qZM 
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where  zu  ■  k  L  .  Reintroducing  the  original  variables  In  accordance 

n  p 

with  equation  (5),  we  have 


T 


(45) 


for  qL  <<  e  .  The  function  j(^)  Is  obtained  by  carrying  out  the 
integral  of  equation  2.60  (Ref.  2)  numerically.  A  second  numerical 
integration  of  equation  (44)  is  then  required  to  obtain  <6>  when  ^ 
not  small.  We  write  the  result  as 


<5y> 

Y 


q(k8L)2 

4k 

w 


cA 

e 


(46) 


where  G  is  plotted  In  Figure  2.  We  note  that  G  decreases  approximately 
linearly  from  one  to  .22  as  Its  argument  varies  between  zero  and  one.  As 
noted  before,  regarded  as  a  function  of  L  ,  <6>  reaches  Its  maximum 
value  at  qL  »  e  ,  yielding 


w  S 

Equation  (46)  can  also  be  conveniently  written  In  the  form 


(47) 


w 


(48) 


The  product  G  O^—)  Is  also  plotted  In  Figure  2.  It  Is  seen  to  have 
maximum  value  of  .27  at  *  .65  •  Thus  we  may  write 
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Extraction  Saturation  Function  for  the  High  Intensity  Regime 


(Cor  Che  optimum  q  ) 


(48) 


The  horizontal  and  vertical  emlttance  requirements  which  are 
applicable  to  this  regime  differ  from  those  derived  In  III.  In  the  case  of 
vertical  emlttance  we  again  have 


but  here  we  require  that  Av^  be  sufficiently  small  to  permit  trapping, 
i.e.,  as  shown  In  (Ref.  2) 


1  + 


Av_  < 


2 


«Y„ 


1  +  a. 


w  Q 


so  that 

2(1  +  a2)  k 

E,  <  - 2 - —  ,/Tji  (50, 

Y  k  8 

w 

In  the  case  of  horizontal  emlttance  It  Is  necessary  to  take  the  coupling  of 
the  transverse  motion  Co  the  optical  synchroton  oscillations  Into 
account.  When  x(o)  does  not  vanish,  equation  (40)  must  be  written 


♦  ■  -e  sin  ♦  +  “  ♦  +  tq*  ~  x(o)]  (51) 

The  condition  for  the  existence  of  a  trapping  potential  at  z  ■  0  Is 
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'  «  *- 


x(o)  <  e  .  Furthermore,  |(o)  -  q  -  (x  -4)  •  -(x  -4)  .  Trapping  for 

x(o)  -  0  requires  (x  -  4)  <  2  JZ  .  In  order  to  take  the  two  effects 
Into  account  simultaneously  we  estimate  (for  a  properly  matched  bean) 


(x  -  4)  x(o)  <  cJ 
o 


Using  equation  (S)  to  regain  dimensional  variables  and  equation  (4)  to 
eliminate  n  •  we  obtain 


•  [(x  4)q  *(o) ^dimensional 


ka(l  +  a£) 


w'  3/2  3/2 

_  £  ■  _ _ _  £ 


4k  y 
w  1 


2k  k 
w  s 


In  conclusion  tie  note  that  these  crude  approximations  may  be 

expected  to  hold  when  QL  »  1  but  jL-  «  I  and  4  «  1  also.  As  seen 

w 

In  Section  IV,  the  principal  effect  of  the  non-linear  terms  Is  to  shift  the 
betatron  frequency  as  the  betatron  amplitude  Increases.  For  the  regime 
discussed  In  this  section,  the  resonance  condition  q  ■  k  plays  no  role, 
the  excitation  of  betatron  oscillations  Is  small,  and  we  expect  the  non¬ 
linear  terms  to  be  less  Important  than  other  omissions  In  the  analysis. 
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VI.  Summary  and  Discussion 

In  Che  preceding  sectiona  we  have  identified  three  distinct 
operating  regimes  for  the  transverse  gradient  wiggler  (TCW) ,  which  we 
characterize  as  low  intensity,  medium  intensity,  and  high  intensity 
regimes.  The  parameters  and  properties  of  these  regimes  are  summarized  in 
Table  I.  In  the  interest  of  simplicity  and  in  recognition  of  the  lack  of 
numerical  precision  in  some  of  our  arguments  we  have  simplified  numerical 
coefficients.  The  more  detailed  description  is  given  in  the  previous  text. 

The  low  and  medium  intensity  regimes  may  be  thought  of  as 
stimulated  Raman  effect  regimes,  in  which  a  quantum  of  betatron  oscillation 
and  an  optical  quantum  are  simultaneously  emitted,  with  the  energy  being 
supplied  by  the  reduction  of  the  amplitude  in  the  transverse  flutter 
motion.  In  more  classical  terms,  the  electron  in  its  longitudinal  rest 
frame  sees  a  ponderomotlve  potential  which  oscillates  at  the  transverse 
betatron  frequency.  This  induces  oscillations  in  y  which,  on  account  of 
the  transverse  field  gradients,  results  in  a  transverse  driving  force  at 
the  betatron  frequency.  The  combination  of  the  flutter  motion  and  induced 
betatron  motion  has  sum  and  difference  frequency  components,  and  as  is 
typical  of  Raman  processes,  only  the  difference  frequency  component  leads 
to  amplification.  The  frequency  of  the  amplified  wave  in  the  laboratory 
frame  is  then  found  by  Doppler  shifting  the  frequency  determined  in  the 
electron's  longitudinal  rest  frame.  As  the  electron  loses  energy  to  the 
optical  wave,  the  center  of  betatron  oscillation  shifts  transversely,  so 
that  the  electron's  average  tranverse  position  shifts  to  a  position  of 


29 


1 


TABLE  I 


Low  Intensity  Regime 

Medina  Intensity  Regime 

High  Intensity  Regime 

definition 
of  regime 

ft 
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« f « ■ 

'  w'  0 

—  »  1 
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16k  .2  L 
w  kg 

k  4 

8  ft  .2 

16k  .2  L 

“k8 

k0  .  .2 

4k  qk0L 
w 

y  Y  /max 

ftf( 

a  \4/3 

v 

k“  ia  i 

18qkw  lv 

A* 

15k 

w 

Lsat 

kl/3 

/  W 

V*' 

% 

a2 

oi 

qke 

E  < 

-v4'3 

4n 

V 

2k  kl2k  1/3 
s  6  w 

4k,k's  2 

ksk8 

E  < 

k8 

k8 

ft3 

h  ' 

2k  k 

8  W 

2k  k 
s  w 

2k  k2k 
s  0  w 

Characteristics  of  the  Various  Operating  Regimes  of  the  Transverse  Gradient  Wlggler. 


*  maximized  with  respect  to  q  for  fixed  L 


weaker  wiggler  field*  The  flutter  motion  is  thereby  reduced  in  amplitude, 
the  reduction  in  kinetic  energy  having  gone  to  supply  that  delivered  to  the 
optical  field.  The  two  regimes  differ  only  in  their  saturation 
mechanism.  Due  to  the  non-linear  nature  of  the  ponderomotive  potential, 
the  amplitude  of  the  driving  force  at  the  betatron  frequency  decreases  as 
the  amplitude  of  the  phase  oscillation  increases,  and  eventually 
vanishes.  Thus  the  driving  force  is  decoupled.  At  the  same  time,  the 
underlined  non-linear  terms  in  the  equations  of  motion  cause  both  the 
betatron  frequency  and  the  driving  frequency  to  shift  as  the  amplitude  of 
the  oscillation  grows.  In  the  medium  intensity  regime  the  intensity  is 

sufficiently  high  that  the  decoupling  effect  occurs  before  the  detuning 

effect  becomes  effective.  In  the  low  intensity  regime  the  growth  rate  is 

slower,  and  the  fact  that  the  driving  force  has  become  non  resonant  becomes 

the  factor  which  limits  growth  of  the  wave.  It  should  be  emphasized  that 
both  regimes  are  non-linear.  The  non-linear  character  of  the  pondermotive 
potential  is  dominant  in  the  medium  intensity  regime,  while  the  underlined 
non  linear  terms  are  dominant  in  the  low  intensity  regime. 

The  high  intensity  regime  may  be  thought  of  as  a  trapped  particle 
regime  analogous  to  the  basic  regime  of  the  high  extraction  variable 
parameter  wiggler  (VPW)  without  transverse  gradient.  The  analogy  is 
especially  close  when  the  parameter  variation  consists  merely  of  a 
reduction  of  the  wiggler  field  intensity  with  z  .  In  both  regimes  the 
particles  are  trapped  in  ponderomotive  potential  wells,  where  they  execute 
optical  phase  oscillations  at  the  optical  synchrotron  frequency.  The 
gradient  in  wiggler  field  shifts  the  phase  center  of  these  oscillations  to 
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one  which,  on  the  average,  causes  a  transfer  of  energy  to  the  optical 
field.  In  the  T GW  case  the  electron  simultaneously  drifts  sidewise  to  a 
region  a  weaker  wlggler  field,  while  in  the  VPW  case  it  convects  to  a 
position  of  weaker  field.  In  both  cases  there  is  a  simultaneous  decrease 
in  the  flutter  amplitude,  thus  accounting  for  the  energy  to  the  optical 
wave. 


The  transition  between  the  medium  and  high  Intensity  regimes 
requires  passage  through  a  region  in  which  the  betatron  and  synchrotron 
frequencies  are  equal.  The  numerical  work  of  Madey  and  Eckstein  suggests, 
as  one  might  expect,  rather  complicated  behavior,  and  oscillations  in 
saturated  gain  may  even  occur.  In  this  connection  we  note  that  the 
saturated  gain  formulas  for  the  high  and  medium  intensity  regimes  match 
badly  at  e  *  1  ,  q  *  k  . 

P 


As  contrasted  to  the  uniform  wiggler  without  gradient  all  of  the 
TGW  regimes  would  seem  to  lead  to  an  Improvement  in  extraction. 

Furthermore,  as  emphasized  in  the  original  proposal  of  the  TGW,  the  energy 
spread  requirements  are  much  less  severe.  Extractions  comparable  to  that 
of  the  VPW  are,  however,  probably  attainable  only  in  the  high  intensity 
regime.  Because  of  the  similarity  of  the  physical  character  of  the  regimes 
one  expects  that  the  optical  intensity  requirements  will  be  similarly 
high.  Indeed,  equation  (50)  for  the  TGW  may  be  written 


a 

.27  - T  k  La  <  .135  k  La 

.  2  w  s  w  s 

1  +  a 
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This  may  be  compared  with 


(‘  •  -28 

for  the  VPW.  This  latter  result  Is  based  upon  equation  4.8  of  Ref.  2 
specialized  to  low  n.  Co  make  It  directly  comparable  to  equation  (50), 

D 

with  ip  ■  *68  rad  and  f^  *  .45  as  Indicated  by  Figure  4.2  b  and  4.2. 

(n  -  n.  at  low  n. ) 
el  d 

By  way  of  direct  comparison,  Madey  has  considered  a  T GW  case  In 
which  Aw  -  5  cm  ,  L  -  20  m  ,  B^nomlnal)  -  6  kG  ,  Y  -  354  ,  and  A  »  1W  • 

For  an  optical  flux  of  10^2  watts/cm^  he  obtains  numerically  an  extraction 
of  14. 22.  (Our  semi-analytlc  approach  yields  12. 82).  For  a  transverse 
magnet  VPW  case  with  the  same  parameters,  except  that  Bw  decreases 
linearly  from  6  to  1.5  KG  as  one  passes  through  the  wiggler,  the  computed 
extraction  (including  both  capture  fraction  and  detrapping  effects)  Is  212. 

The  effective  vertical  focussing  forces  In  the  TGW  are  stronger 
than  those  In  the  VPW,  a  fact  which  tends  to  make  the  emittance 
requirements  more  restrictive.  On  the  other  hand,  for  the  VPW  with 
constant  the  trapping  potential  (i.e. ,  $Ym)  *s  reduced  by  the  factor^ 

r( )  •  The  two  effects  approximately  cancel  for  the  comparison  example 

considered  above. 

For  the  VPW  with  transverse  magnet,  the  horizontal  emittance 
restriction  is  very  weak,  as  there  are  negligible  horizontal  focussing 


forces  and  negligible  horizontal  field  variation*  For  convenience  of  beam 
design  one  takes  -  Ey  ,  but  g^  »  g^  would  presumably  be  achievable 
if  the  total  emittance  constraint  were  found  to  excessively  restrictive* 

In  contrast,  for  the  example  being  considered,  the  horizontal  emittance 
requirement  equation  (50)  is  more  than  fifty  times  more  restrictive  than 
the  vertical  one.  Thus  to  take  full  advantage  of  the  vertical  beam 
emittance  available,  a  ribbon-like  beam  structure  would  be  required,  and  it 
may  prove  very  difficult  to  obtain  such  a  small  horizontal  emmltance. 

Although  we  have  not  carried  out  the  parameter  exploration  that 
would  be  required  to  reach  a  decisive  conclusion,  the  above  discussion 
provides  a  strong  indication  that  compared  on  a  steady  state  operation 
basis,  the  VPW  will  prove  substantially  superior  to  the  TGW  for  high 
extraction  single  pass  operation. 

An  important  issue,  which  we  have  not  studied,  is  that  of  the 
build-up  of  oscillations  from  a  low  level,  and  the  nature  of  the  optical 
pulse  form  which  emerges  when  the  system  is  driven  by  the  electron  micro- 
bunches  typically  produced  by  an  RF  LINAC.  Qualitative  changes  in  the 
operating  regimes  take  place  as  build-up  proceeds,  and  both  devices  are 
subject  to  detrapping  by  field  fluctuations.  It  is  entirely  possible  that 
Important  differences  will  emerge,  and  that  either,  neither,  or  both  may 
prove  desireable  from  this  point  of  view. 


In  conclusion  we  wish  to  acknowledge  that  the  development  of  the 
preceding  theory  was  facilitated  by  the  numerical  work  of  Eckstein  and 
Madey.  Their  work  has  also  provided  us  with  the  opportunity  to  make  some 
detailed  quantitative  and  qualitative  comparisons  with  their  numerical 
results.  The  agreement  is  generally  very  good,  and  we  have  found  no 
significant  discrepancies  in  regions  where  we  consider  that  our  theory 
should  apply. 
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